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Abstract 

The article deals with the description of the statistical behavior of Gaussian packets on a metric graph. 
Semiclassical asymptotics of solutions of the Cauchy problem for the Schrodinger equation with initial data 
concentrated in the neighborhood of one point on the edge, generates a classical dynamical system on a graph. 
In a situation where all times for packets to pass over edges ("edge travel times") are linearly independent 
over the rational numbers, a description of the behavior of such systems is related to the number-theoretic 
problem of counting the number of lattice points in an expanding polyhedron. In this paper we show that 
for a finite compact graph packets almost always are distributed evenly. A formula for the leading coefficient 
of the asymptotic behavior of the number of packets with an increasing time is obtained. The article also 
discusses a situation where the times of passage over the edges are not linearly independent over the rationals. 

1 Introduction. 

In this paper we study semiclassical solutions of the Cauchy problem for a time-dependent Schrodinger equation 
on a metric graph. This article is, in a s certain sense, a continuation of [5], [3] where you can find necessary 
definitions and formulations (nevertheless, we quote all necessary terms here, in section ll~2")) . There (and in [3]) 
are also references to some articles and reviews related to the study of differential equations on metric graphs. 

In [3] a formula for the asymptotics of the number of Gaussian packets was written for a star graph, and a 
question of finding the leading coefficient for the case of an arbitrary finite compact graph was raised. In this 
paper that problem is solved and a general formula for that coefficient is obtained (see Theorem l2.2p . Moreover, 
in [3J the uniformity of distribution of packets for the special case of two vertices connected by three edges 
was proven. In this paper we present the proof of uniformity of the asymptotic distribution of packets over an 
arbitrary finite compact graph for almost all edge travel times (see Theorem 12. ip . 

In the second part (section [3J we discuss a description of the asymptotic behavior of packets in the case 
of linearly dependent over Q edge travel times. In this situation there is no correspondence (see [2J, [5J for 
details) between the number of packets and the number of nodes of an integer lattice that lie on some faces of 
an expanding simplex, but analysis of their number is still possible. If rank of the system of edge travel times 
(lengthes) is equal to one then the number of packets grows only in a finite interval of time and stabilizes at a 
certain value, which depends on the lengths of cycles. Also we consider an example of a graph, whose rank of 
the system of lengthes equals two. 
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1.2 Terms and Definitions. 

Recall (see, e.g., [5], [5] and references therein), that a metric graph is a one-dimensional cell complex whose 
edges are parameterized curves. We denote a geometric graph by T, its edges by jj, its vertices by aj. A set 
of all edges adjacent to the vertex a we denote by T(a). We consider only finite metric graphs. Let E and V 
stand for the number of edges and the number of vertices respectively. 



Let Q be an arbitrary real valued continuous function on T, smooth on the edges. Schrddinger operator 



Hi> = -h 2 ^^- + Q(xW(x) (1) 
is defined on the set of functions from Sobolev spaces ip £ ®J2-H 2 (lj)i satisfying the following boundary 

3 

conditions at the vertices: 

1. function ip is continuous on T; 
2. 



E a ^M = °i ( 2 ) 



i 3 ev(a m ) 

in all internal vertices (i.e., the vertices of degree greater than one); 
3. ip{a m ) = in all the external vertices, i.e. vertices of degree one. 

Here otj = 1 for each edge emerging from the vertex, and ay = — 1 for each incoming edge. 

These conditions are called natural (see |10j). They, in particular, ensure self-adjointness of the operator H. 

A time- dependent Schrddinger equation on the graph T is an equation of the form 

ih^ = H1>, (3) 

where a semiclassical parameter h — >■ 0. We choose initial conditions that have the form of a narrow packet 
localized near the point xq, which lies on the edge of the graph: 



Sq(x) = a(x — xq) + b(x — xq) + c 



(4) 



where b and c are real constants, and a and K are complex. The imaginary part of a is positive. Normalization 
factor h^ 1 ' 4 is introduced to ensure that the initial function ip(x,Q) is of order one in the norm of L 2 (T). Due 
to the positivity of the imaginary part of a the initial function is localized in a small neighborhood of Xq- 
ip(x,0) = 0{h°°) with \x — Xq\ > S > (S is independent of h). For simplicity we assume that there are no 
turning points (see, e.g., [7J, [5]) on the V ; their presence can be accounted for in the standard way (see [7]). 

Asymptotic solution of the Cauchy problem ©-(jU) is described in [3], the explicit formulas are given therein. 

Theorem 1.1 (See J^/) Solution of the Cauchy problem\3\ -\J\for t G [0, T] (T does not depend on h), is given 
by the following formula 

N(t) 

i>(x,t,h) = h-^ipjity^^^ + OiVh), (5) 

where the functions <Pj(t), Sj(t) are explicitly expressed in terms of the solutions of two hamiltonian systems. 

Each term in the sum (O is localized in a small neighborhood of Xj{t). Here we assume that all the terms 
that are localized in the same point Xj(t) form one Gaussian packet. Later, under the N(t) we would mean the 
number of such packets. 

In [3j it is shown that in moments of penetration of the vertices of the graph a quantum packet is divided 
into m packets (m is a degree of the vertex): one is reflected and m — 1 are "diffused". The amplitude of the 
packet is divided as (m — 2) : 2 (2 for the scattered and (m — 2) for the reflected packet). 

We consider the asymptotical behavior of function ip(x, th) as t^-oo. Namely, we will see how the number 
of Gaussian packets N(t) changes in time. Note that this problem differs from the task of describing the 
asymptotic solution of the Schrodinger equation at t— >oo, as the error estimation is valid only for finite times. 
From a physical point view, this means that we are considering big t, but much smaller than 1/h. 

Let tj stand for j-th edge travel time (this is an analog of length) . 



Lemma 1.1 Travel time of any edge of the graph depends only on the initial data and is the same for any 
Gaussian packet on each fixed edge. 

Proof. Let us take an edge with index j. By construction a solution on each edge (using the method of 
Maslov complex germ), we find that P 2 + Q(x) = Ej holds (by the energy conservation law for Hamiltonian 
systems), where the value of Ej is determined by the initial conditions. Writing P as ^r, we obtain for edge 

b 

travel time on the edge an explicit expression tj — J — ^===, where the integral is taken over the edge. The 

value of Ej is the same for all edges, as the potential Q(x) is assumed to be continuous function in the vertices, 
and P can only change sign (as demonstrated by a construction of the solution). So the formula Ej = P 2 + Q(x) 

determines the same value Ej = E for all packets. At the initial moment of time E — (^f) 2 + Q(xo)- Values 
of tj may be different, since the restriction of the potential Q(x) on the edges may be different. 

Definition. Let us consider the number of packets coming out of a fixed vertex to a fixed edge. The leading 
coefficient of the asymptotics of this number is called a radiation coefficient. 

Correctness of this definition will be shown in the proof of the theorem 12.11 



2 The case of linearly independent over Q edge travel times. 
2.1 Main results. 

Theorem 2.1 (About uniformity of distribution) Consider a finite connected graph T. Suppose that for 
any vertex its degree is not equal to two. Suppose that there are no turning points on edges. Consider an edge 
e. Let F e (t) be a ratio of N e (t) (number of packets on the edge e) to N(t) (total number of packets). Then for 
almost all incommensurable (i.e. linearly independent over Q) numbers t\. . . .tg, a ratio of F e (t) to the length 

( I \ : 

of e tends to a constant }, t,- I as time increases. 

\h V 

Remark 2.1 It means that the distribution of number of packets tends to a uniform distribution as time in- 
creases. 

Proof. Let us choose on any edge with travel time tj a segment dg with travel time r. Let us find N T (t)/N(t) . 
We know (see [5]) that N(t) = Ct E ~ x + o(t B ~ 1 ). Let us find N r (t). Since the number of packets changes only 
in vertices and there are no turning points, then: 

N T (t) = #_>«,(*) - ~t) + N^ g (t) - N^ g (t - t). (6) 

Here N^d{t) stands for the number of packets which arrived at the edge from point d. 

It is clear that the number of packets arrived at point d at time t equals the number of packets, which came 
out from the nearest vertex a at time t — T\. Here T\ is a travel time from a to d. By N a ^d(t) we denote the 
number of packets, which came from a to d. 

We have to know asymptotics of the number of packets that come out of a vertex a. Packets can come out 
of the vertex only at times that are linear combinations (with nonnegative integer coefficients) of edge travel 
times. 

The number of release moments (when at least one packet comes out of the vertex a) is described by the 
number of set {nj} satisfying inequations of a kind: 

nifji + . . . + n m ti m <t, (7) 

where tj is a travel time of the j-th edge. 

Since the leading part of asymptotics of the number of packets is defined by the volume of a simplex defined 
by ([7]), events with maximal numbers of summands happen more often. In other words, packets arrived at our 
vertex should visit all edges. I.e. 

N a ^ d (t) =R a t E +o(t E ). (8) 

For almost all t\,...tE the estimation can be improved (see [T]). There exists K a such that N a ^d(t) = 
R a t E + K a t E ~ 1 +o(t E ~ 1 ). Let us show that R a , which is called a radiation coefficient, does not depend on the 



choice of a vertex. Consider vertices a and b. There exists a path connecting a and b. Let 5 be its travel time. 
Any packet coming out from a to d over time that does not exceed 26 generates at least one packet that come 
out from b to d! . This is correct for packets coming out from b. We obtain inequations: N a ^ d {t + 26)>Ni } ^ ) .d' (t) 
and Nb-+d'(t + 28)>N a ^ d (t). We know that N a ^ d {t) = R a t E + o(T E ), N b ^ d ,(t) = R b t E + o{T E ). Thus 
R a t E + o(t E ) = R b t E . Hence R a = R b . 
Let us modify the expression for N T (t) 

N T (t) = R(t ~ T X ) B + K a (t - Ti)^- 1 - R(t - Ti - t) e -K a (t-T x - r)^ 1 + R(t - T 2 ) E + K b {t - T 2 ) E - 1 - 
R(t -T 2 - t) e - K b (t - Ti - t)^- 1 + o^- 1 ) = 2ERrt E - 1 + o^- 1 ). 
Thus we obtain 

MA (9) 
N(t) C () 

It remains to show that a coefficient near r has the required form. 

We consequentially take edges as dg and then summarize obtained expressions: 

Hence, 

E 

C = 2ER^2tj- (11) 

The proof is completed. 

Remark 2.2 For any two vertices numbers K a and K b are related with inequation \K b — K a \<ESR. Here 8 
is a travel time for any path from a to b. 

In the proof we obtain the following statement: 

Consequense 2.1 (Relation between coefficients C and R) The leading coefficient for the number of pack- 
ets C and the radiation coefficient R, for almost all edge travel times, are related in the following manner: 

E 

C = 2ERj2 t 3- (12) 
i=i 

Theorem 2.2 (About the leading coefficient of the number of packets) Consider a finite connected com- 
pact graph T. Suppose that there are no vertices of degree two. Suppose that there are no turning points on 
edges. Then for almost all incommensurable numbers t\ , . . . the leading coefficient has the following form: 

E 

i S* j 

Proof is based on (fT3| and the following lemma. 

Lemma 2.1 Let us consider a finite connected graph with incommensurable edge travel times U(i = 1 . . . E) 
and a number of independent cycles f3. Let B be an arbitrary vertex. Then for almost all edge travel times the 
number of packets arriving at B at time T asymptotically equals 

R(T) — T E . 

2 E E\Y[t J 



Remark 2.3 This is equivalent to 

R 1 1 



2 V ~ 1 E\ e ■ 

3=1 



Proof. Is is sufficient to consider paths of packets that traveled upon ali edges. Oniy those paths give us 
the ieading coefficient. Let A be an initial vertex. For each such path we can construct a "code" i.e. a sequence 
of coefficients of a corresponding chain with coefficients in Z2 . It is ciear that the code does not change under 
a path homotopy. Let us find the number of all possible codes. Consider cross connections, i.e. edges that are 
not in the spanning tree. Parity of passages on cross connections defines a path's homotopy class. All chain 
coefficients are defined by coefficients on cross connections. Thus the number of possible codes equals 2^. Now 
for every code (ci, . . . , ce), Cj € {0, 1} we associate times 



J2u(c, + 2m)\n % eNU{0} 



At every such time (they are different) at least one packet arrives at the vertex B. The number of such times 
that are less than T asymptotically equals to 



2 E Eltf-t E ' 



Finally we summarize this over all possible codes. The proof of the lemma is finished. At the end we apply 
Euler's relation (3 = E — V + 1 (see, for example, [S]). 



3 The case of linearly dependent travel times. 

In this section we assume that travel times are linearly dependent over Q. It means that there is no one-to-one 
correspondence (described in [2], [3]) between the number of packets and the number of integer lattice points 
in an expanding simplex. 

For the simplest example consider a star graph with three edges of the same length. The number of packets 
reaches three and does not increase. While the number of integer lattice points grows with time. 

Statement 3.1 Consider a finite graph with edge travel times t\,= nito . . . ,t E = nEto, where m S N and 
GCD(m,...,n E ) = 1- 

Then starting from a certain time, the number of packets becomes constant. 
1) If there exists a cycle with travel time 2kto, k € N then 

E 

N(T)=2Y / n l . 

2) Otherwise 

E 

N{T) = Y,n i . 
i=i 

We omit the proof. 

Statement 3.2 Consider a star graph with three edges &i,&2,e§ with travel times t\ = ntQ,t2 = mto,t3, where 
n G N, m G N (GCD(n.m) — I), and £3 is such that rank {^1,^2,^3} over Q equals 2. Then the number of 
packets asymptotically equals 

._. T ( m + n 1 \ 



Proof. The number of packets changes only at times of a kind 2((na + mf3)to + £37), where a, (3, 7. Note 
that this representation is not unique. If we can present a time moment with a > 1, (3 > l,j > 1, then at that 
moment N(T) does not change (since packets arrive at a vertex from all three edges) 

Remark 3.1 Number na can be represented as na' + m/3' fl3' ^ 0) if and only if a> m. 

Hence to define asymptotics of N(T) it is sufficient to consider the following times: 

1) j ^ 0, /3 ~ 0,1 < a < m. At those times N(T) increases by 1. The number of such times equals the number 
of pairs (a, 7) such that 2(ato + jts) < T. The number asymptotically equals: (m — 1)^-- 

2) 7 7^ 0, a = 0, 1 < (3 < n. Similarly the number of pairs asymptotically equals (n — 1)^-- 

3) a = 0,/3 = 0. The number of such times that are less than T increases as Here the number of packets 
increases by 2. 

4) 7 — 0. Let us show that linear combinations an + f3m contain all natural numbers that are greater than a 
certain fixed number. 

Lemma 3.1 If (to, n) = 1 then there exists M such that any N > M can be presented as N = an + /3m, where 
(3>l,a> 1. 

Proof. Numbers to, 2to, . . . , nm give all residues modulo n: 

{(3m}'^ =1 = {fam = i + nk t \i e [0,n- l] 1 k l <= Z+}. 

Let fco = max fcj. Choose any N > nko. Then N = i + nk (k > ko, < i < n) and N ~ ftm = n(k — fcj). Thus 
any number that is greater than nko can be represented as an + (3m, where 1 < (3 < n,a > 1. 

Asymptotics of the number of such times coincides with asymptotics of numbers 2kt n (k e N) that are less 
than T. Since at such times 2 packets arrive from e\ and e 2 , then N(T) increases by 1. 
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